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We find that the overlapping of a topological quantum color code state, representing a quantum 
memory, with a factorized state of qubits can be written as the partition function of a 3-body 
classical Ising model on triangular or Union Jack lattices. This mapping allows us to test that 
different computational capabilities of color codes correspond to qualitatively different universality 
classes of their associated classical spin models. By generalizing these statistical mechanical models 
for arbitrary inhomogeneous and complex couplings, it is possible to study a measurement-based 
quantum computation with a color code state and we find that their classical simulatability remains 
an open problem. We complement the meaurement-based computation with the construction of a 
cluster state that yields the topological color code and this also gives the possibility to represent 
statistical models with external magnetic fields. 



PACS numbers: 03.67.Lx, 03.67.-a, 75.10.Hk, 05.50.+q 



I. INTRODUCTION 

Recently, a very fruitful relationship has been estab- 
lished between partition functions of classical spin models 
and a certain class of quantum stabilizer states with topo- 
logical protection [H , • The topological quantum code 
states considered so far in these studies correspond to the 
toric code states introduced by Kitaev 0], [J. The clas- 
sical spin model that emerges when a planar toric code is 
projected onto a product state of singlc-qubits with very 
specific coefficients is the standard classical Ising model 
in two dimensions with homogeneous real couplings and 
zero magnetic field. 

Single-qubit measurements also appear naturally in a 
measurement-based computation (MQC) scheme 0- 
Thus, these connections between classical spin models 
and topological quantum states are also useful to test 
whether those topological states are efficiently classically 
simulable with MQC. It has been shown that MQC with 
a planar Kitaev code state as input can be efficiently 
simulated in a classical computer if at each step of the 
computation, the sets of measured qubits form simply 
connected subsets of the two-dimensional lattice [1] . The 
connection of classical spin models with measurement- 
based quantum computation has been shown to be useful 
to prove the completeness of the classical 2D Ising model 
with suitably tuned complex nearest-neighbor couplings 
in order to represent the partition function of the classical 
Ising model on arbitrary lattices, with inhomogeneous 
pairwise interactions and local magnetic fields [7[. 

Topological color codes (TCC) were introduced to im- 
plement the set of quantum unitary gates of the whole 
Clifford group by means of a topological a stabilizer code 
in a two dimensional lattice [8|, and then generalized 
to three dimensional lattices in order to achieve a uni- 
versal set of topological quantum gates Q. These 2D 
and 3D realizations of TCC are instances of general D- 
dimensional realizations. We call those lattices related 
to these codes as D-colexes (for color complexes), and 
they are D-dimensional lattices with coordination num- 
ber D + 1 and certain colorability properties. Moreover, 



this codes can also appear as the ground state of suitable 
Hamiltonians, and the corresponding quantum systems 
are brane-nct condensates. [1C| . 

Given these nice properties exhibited by the topolog- 
ical color codes, it is natural to ask what type of clas- 
sical spin models can be constructed out of them and 
see whether they belong or not to the same universality 
class of the classical Ising model arising in the Kitaev 
model. In this work we address this issue and find that 
the the overlapping of a TCC state with a product state 
of single qubits with appropriate coefficients is mapped 
onto the partition function of the 3-body classical Ising 
model on the dual lattice of the original lattice where the 
color code is defined. For concretcness, we consider the 
triangular and the 'Union Jack' lattice for these classical 
many-body spin systems. This represents a sharp differ- 
ence with the result obtained with the topological states 
in the Kitaev code. In fact, the universality classes of the 
3-body classical Ising model in several lattices are quite 
different from the corresponding universality class of the 
standard 2-body Ising model. 

Moreover, we also study the topological color code 
states in a MQC scenario to test their classical simu- 
lability. We find that the current state of knowledge in 
statistical mechanical models with 3-body interactions, 
arbitrary inhomogeneous complex couplings and lattice 
shapes is much less developed than the 2-body Ising 
model which is relevant for the case of the toric code 
states. Thus, we conclude that the classical simulability 
of TCC states with MQC remains an open problem. 

In a MQC, the usual initial many-particle entangled 
state is a cluster state Q , instead of a topological code. 
Then, we also show how to construct a color code state 
from a certain cluster state. Interestingly enough, this 
construction turns out to be useful for the description of 
statistical mechanical systems with 3-body interactions 
in the presence of an external field. 

This paper is organized as follows: in Sect|II]we give an 
introduction to the topological color code states needed 
to present in Sect lllll the mapping onto the classical 3- 
body Ising model in the triangular and Union Jack lat- 



Figure 1: An example of 2-colex. Both edges and face are 3- 
colorable and they are colored in such a way that green edges 
connect green faces and so on and so forth for red and blue 
edges and faces. 

tices. In Sect llVl we study the measurement-based quan- 
tum computation with topologocial color code states by 
generalizing the results of the previous section. In SectlVl 
we show how to prepare a topological color code from a 
cluster state as those introduced in MQC. This is also 
useful for studying partition functions of statistical me- 
chanical models with 3-body interactions and external 
magnetic fields. ScctlVIlis devoted to conclusions. 



II. TOPOLOGICAL COLOR CODES 

A. Construction 

Let us start by recalling the notion of a Topological 
Color Code (TCC) in order to sec what type of classi- 
cal spin models we obtain from them with appropriate 
projections onto factorizcd quantum states and specific 
lattices. 

A TCC, denoted by C, is a quantum stabilizer error 
correction code constructed with certain class of two- 
dimensional lattices called 2-colexes. The word colex is 
a contraction that stands for color complex, where com- 
plex is the mathematical terminology for a rather general 
lattice [ll|- A 2-colex, denoted by £2, is a 2D trivalent 
lattice which has 3-colorable faces and is embedded in 
a compact surface of arbitrary topology like a torus of 
genus g. A trivalent lattice is one for which three edges 
meet at every vertex. The property of being 3-colorable 
means that the faces (or plaquettes) of the lattice can be 
colored with these colors in such a way that neighbor- 
ing faces never have the same color. We select as colors 
red (r), green (g) and blue (b). An example of a 2-colex 
construction is shown in Fig. [T] 

Edges can be colored in according to the coloring of 
the faces. In particular, we attach red color to the edges 



Figure 2: An hexagonal lattice is an instance of 2-colex. Num- 
bered vertices belong to the face f. Vertices labeled with let- 
ters correspond to the red string 7 displayed. 7 is an open 
string, because it has an endpoint in a red face. 

that connect red faces, and so on and so forth for the blue 
and grcc edges/faces. When studying higher dimensional 
colcxes, it turns out that the coloring of the edges is the 
key property of D-colexcs: all the information about a 
ZJ-colex is encoded in its 1-skeleton, i.e., the set of edges 
with its coloring j7(]|. 

Given a 2-colex £2, a TCC C is constructed by placing 
one qubit at each vertice of the colored lattice. Let us 
denote by 2J, € and ^ the sets of vertices 0, edges e 
and faces f, respectively, of the given 2-colex. Then, the 
generators of the stabilizer group, denoted by 5, are given 
by face operators only. For each face f, they come into 
two types depending whether they are constructed with 
Pauh operators of X- or Z-type: 

(1) 

:=(g)^o, 

and there are no generators associated to lattice vertices. 
For example, an hexagonal lattice is an instance of a 2- 
colcx, see Fig. [21 The operators for the face f displayed 
in the figure take the form Aj = A1A2A3A4A5A6, = 
Z1Z2Z3Z4Z5Z6. A given state l^^c) e C is left trivially 
invariant under the action of the face operators, 

Af|*c) = |*c), ZflvPc) - |*c), Vfe^. (2) 

An erroneous state |'I')o is one that violates conditions 
([2]) for some set of face operators of either type. As the 
generator operators Xj, Zf G S satisfy that they square to 
the identity operator, (Aj)^ = I = (Zf)^,Vf £ ^, then an 
erroneous state is detected by having a negative eigen- 
value with respect to some set of stabilizer generators: 

Af|«')e = -|*)e and/or Zf|*)e = -|*)e. 

Interestingly enough, it is possible to construct a quan- 
tum lattice Hamiltonian such that its ground state is 
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degenerate and corresponds to the TCC C, while the er- 
roneous states are given by the spectrum of excitations 
of the Hamiltonian [3| . Such Hamiltonian is constructed 
out of the generators of the topological stabilizer group 
S, 

H, = -J2iXf + Zf). (3) 

feS 

The ground state of this Hamiltonian exhibits what is 
called a topological order [l^l , as opposed to a more stan- 
dard order based on an spontaneous symmetry breaking 
mechanism. One of the signatures of that topological or- 
der is precisely the topological origin of the ground state 
degeneracy: the number of degenerate ground states de- 
pends on topological invariants like Bctti numbers p^ . 
In two dimensional lattices, the relevant Betti number 
corresponds to the Euler characteristic x of the surface 
where the 2-colex is embedded. 



B. String- net operators 




Figure 3: A 4-8 lattice is an instance of 2-colex. A closed 
string-net is displayed, composed of 3 strings of different col- 
ors meeting at a branching point. The string-net is closed 
because at each face we find an even number of its vertices. 



In order to better understand both the ground state 
and excitations of this Hamiltonian and their topolog- 
ical properties, it is rather convenient to introduce the 
set of string operators that can be defined on a 2-colex 
£2- String operators are generalizations of face operators 
(HI) that can be either open or closed, i.e., with or with- 
out end-points. These strings are topological and like 
in Kitaev model, the homology is defined on Z2 since 
we work with two-level quantum systems located at the 
sites of the lattice. However, in a TCC we have an addi- 
tional ingredient to play around: color. Let us split the 
sots of edges and faces into colored subsets denoted by 
IS := (£r U (£g U (Bb and ^ := dr^Sg^Sb, where €r is the 
subset of red edges, and similarly for the rest of subsets. 

A colored string 7 is a collection of edges of a given 
color. Thus, a blue string 7 takes the form 7 = {ci} 
with Ci £ The definition of colored string operators 
is completely analogous to that of face operators: 

X^-.^^X,, Z^:=(g)Ze, (4) 

where, in turn, X^ = X„-^ ®X„2 if fi and O2 are the sites 
at the ends of the edge e, and similarly for = Z„-^ iS) 
Z„2- For instance, consider the red string operator in 
Fig. [21 where we have X^ = XaX^XcX^XeX f ■ ■ ■ , Z^ = 

ZaZbZcZdZ^Zf ■ ■ ■ . 

Colored strings are open if they have endpoints. These 
endpoints are localized at faces which share color with 
the string. In particular, a face f is an endpoint of 7 if 
the number of edges of 7 meeting at f is odd, see Fig. [21 
In terms of string operators, a face f is an endpoint of 7 
if {Xj,Zf} = or, equivalently, if {Z^,Xf} = 0. Thus, 
open string operators do not commute with those face 
operators in their ends. In other words, a string op- 
erator that commutes with all the face operators must 



correspond to a closed string, that is, a string without 
endpoints. In terms of the Hamiltonian Q, string op- 
erators produce quasi-particle excitations at their ends 
when applied to the ground state. These quasi-particle 
excitations are Abelian anyons. 

Closed strings are mainly of two types. They can be 
homologically trivial, meaning that they are the bound- 
ary of certain area of the surface, or homologically non- 
trivial. In terms of operators, a string is a boundary iff 
its string operators belong to the stabilizer group S of 
the color code C. Such boundary string operators are 
thus products of face operators. In fact, face operators 
themselves are the basic boundary string operators. 

Although we have introduced colored strings and the 
corresponding operators for illustrative purposes, in fact 
in a TCC we have to deal with more general objects, 
namely string-nets. A string-net is a collection of strings 
meeting at certain branching or ramification points. An 
example of these types of configurations are shown in 
Fig. [3 

A string-net 7 is a collection of vertices 7 C 5J. Equiv- 
alently, 7 is a formal sum of lattice vertices S QJ with 
coefficients 7,, G Z2, i.e., 

7 = XI 

ve'S 

where 70 = 1 if G 7 and 7d = otherwise. Given a 
string-net 7, we define the string-net operators 

X^:=(g)X'^", Z^:=(^Z'" (6) 

D 

Just as in the case of colored strings, we can talk about 
open and closed string-nets, and about trivial and non- 
trivial closed string-nets. In terms of operators, the sit- 
uation is exactly the same as with strings. That is, a 
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string-net 7 has an endpoint at a face f if {Xj, Zj} = 0, 
it is closed if its string-net operators commute with aU 
the face operators, and it is a boundary if it is a prod- 
uct of face operators. In order to translate these ideas 
into purely geometric terms, we can define the boundary 
operator 



del := X]^ff' 



7 n f I is even, 
7 n f I is odd, 



(7) 



where |7nf| is the number of vertices that 7 and f share. 
Thus del is the formal sum of the endpoints of 7. It is 
also natural to define an operator dc for faces 



dcf := ^.T„D, 



0, O^f, 

1, D e f. 



(8) 



so that 9cf is the string-net composed of the vertices of 
f. With this definitions, 7 is closed if and only if 



del = 0, 



(9) 



and it is a boundary if and only if there exist a collection 
of faces S = 5ff such that 



7 = deS. 



(10) 



It is possible to give explicit expressions for the states 
of the TCC or, equivalently, for the ground states of the 
Hamiltonian ^ . The states are superpositions of all pos- 
sible closed string-nets, a typical feature of t he g round 
states of systems with topological order [l^, [iSj The 
following is an un-normalized ground state for any given 
2-colex H, [l3l 



«>|aj| 



\^e) :=n(l + ^t) |0) 



= 5: X, |0)«l^l =: Y: l7>, 

76ro 76ro 



(11) 



where |QJ| is the number of vertices in the 2-colex £2, 
Tq denotes the set of boundary string-nets and |0} is the 
eigenstate Z|0) = |0). 

The degeneracy of the ground state or, equivalently, 
the number of logical states encoded in the color code, de- 
pends on the topology of the lattice. For a general 2-colex 
£2 with Euler characteristic x(22) |2J| - |£| + the 
number k of encoded qubits is given by fc = 4 — 2x(£2) := 
2hi Q, where /ii is the first Betti number of the surface 
where the 2-colex is embedded [l3l. These additional 
ground states can be obtained from the one given by (jlip 
by the action of the encoded logical operators Xi , Zi with 
I = 1, . . . , fc. These, in turn, take the form of string-net 
operators of non-trivial closed string-nets 0. 

For the purpose of this work, we shall be interested 
only in a representative ground state like (jlip . Thus, we 
will have to consider suitable surface topologies such that 
the corresponding TCC is unique. We will return over 
this issue later when we consider particular lattices. 



III. CONNECTION WITH CLASSICAL SPIN 
SYSTEMS 

A. Overlap and Partition Function 

Now, we come to the issue of what type of classical spin 
models may arise from the color code state (fTT|) when we 
project it onto a product state of a number of qubits given 
by |QJ|. In this section we shall not consider the most 
general factorizcd state, but one specifically adapted for 
the purpose of this connection in its most simple form, 
namely. 



l^p) := <S> l'^)"! I-^)" ■= cosh(/3J)|0)„ + sinh(/3J)|l)„, 

(12) 

with f3 := l/k^T the inverse temperature parameter. 

The classical spin model arises when computing the 
overlapping between the ground state of the color code 
Hamiltonian pT|) and this factorized state (fT2|) . 



0{(3J) := (*e|$p). 



(13) 



Using pT|) and we get the following expression for 
this overlapping. 



7ero 



(7l(X)l0)o = (cosh(/3J)) 



m ^ ,i7i 

76ro 



(14) 

u :— tanh(/3J) and I7I is the number of vertices of 7. 

We want to relate to the partition function of a 
classical spin system. So let £2 be an arbitrary 2-colex. 
Consider the dual lattice A. The vertices of A correspond 
to the faces of £2, and the faces of A are vertices in 
£2- In particular, A is a lattice in which all faces are 
triangular and vertices are 3-colorable. Moreover, for any 
such lattice A there exist a suitable dual 2-colcx £2. 

So let us associate a classical system to A by attach- 
ing classical spin variables ai = ±1 to each of its sites i 
(equivalently, to each face f of £2). The classical Hamil- 
tonian is 



H := —J ^ (JiajCTk, 



(15) 



where J is a couphng constant and the sum X](ijfe) 
is over all triangles with spins aiGjOk at their vertices. 
Thus, we have a classical Ising model with 3-body inter- 
actions. The case J > corresponds to a ferromagnetic 
model with an even parity to be discussed below, and 
similarly J < to an antiferromagnetic model with odd 
parity. The partition function of the model is 



Z{PJ) 



(16) 



where the sum X]{cr} over all possible configurations of 
spins. The point then is that we have 



Z(/3J) = 2^^0(/3J), 



(17) 



Figure 5: A typical chain of triangles S £ Aq in a triangular 
lattice. It is understood that only part of the lattice is dis- 
played. Black triangles represent the elements of S. The fact 
that 6 € Ao means that at each vertex always meet an even 
number of triangles. 



Figure 4; Two instances of dual lattices of a 2-colex, which 
have triangles as faces and have 3-colorable sites. The trian- 
gular lattice (a) is dual to the hexagonal one. The Union Jack 
lattice (b) is dual to the square-octogonal or 4-8 lattice. 



where <5 = j fc) ^ijk^ijk is a chain of triangles, that is, 
a formal sum over triangles with binary coefficients, and 
\5\ is the number of triangles in 5. Using the identities 



where N is the number of sites. 

Before we show why this identity holds, let us give a 
pair of representative examples of dual lattices £2 and 
A. First, if the 2-colex is an hexagonal lattice then the 
dual lattice A is a triangular lattice, Fig. |4^a). Second, 
if the 2-colex is a square-octogonal lattice (also denoted 
by 4-8 lattice), then its dual is a Union Jack lattice, see 
Fig. mb). The relevance of these examples is two- fold. 
On the one hand, the hexagonal lattice is the simplest 
lattice for a 2-colex and the 4-8 lattice is the simplest one 
when we want to obtain TCC with certain transversality 
properties for quantum computation (see below). On the 
other hand, 3-body classical Ising-models on both lattices 
have been studied in statistical mechanics to some extent. 

To prove (fT7|) , let us start expanding the partition func- 
tion Z{(3J) using the following identity, 

^(3Ja,a,a^ = cosh(/3J) -f- a,(jjak sinh(/3J). (18) 

Inserting it into (|16p , we may expand the partition func- 
tion as 

Z{PJ) ^ (cosh(/3J))^ E n (1 + ^'^■^'^j'^k). (19) 

Let us rewrite (|45|) in the form 

Z(/3J) = (cosh(/3J))^^Hl^l5] n ('^^^J-^fc)''^"' 

(20) 



= 0, ^ a"^ ^ 2, (21) 

f7=±l CT = ±1 

where Uq and rig are odd and even numbers, respectively, 
we get 

Z{f3J) = (2cosh(/3J))^ (22) 

(SgAo 

where Ag contains those chains of triangles such that at 
any given site i an even number of triangles meet, as 
shown in Fig. [5l In fact, this type of expansion is called 
a high-temperature expansion of the partition function 
of a statistical mechanical model [l^ . 

In order to compare (j40p and (j52[) . we simply observe 
that triangles A^-fe in A correspond to vertices of the 2- 
colex G QJ. This correspondence relates in an obvious 
way a string-net 7 with a triangle chain 6, in such a 
way that Ag is identified with Tq. Therefore, we have 
the desired relationship between the overlapping and the 
partition function p7p . 

Although the previous derivation was performed for a 
model with uniform couplings, it is possible to obtain a 
completely analogous result for triangle dependent cou- 
plings Jijk- For simplicity we have preferred to do the 
exposition with uniform couplings because, in fact, the 
case of non-uniform couplings is contained in the more 
general case of non-uniform couplings with non-uniform 
external field, to be considered in section IVl 
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B. Consequences 

We hereby draw a series of very important conse- 
quences from these resuhs, that will continue in the next 
section. 

i/ Interactions: We see that there is a clear qualita- 
tive difference between topological color code states and 
Kitaev's toric code states since they yield quite different 
type of spin interactions: a TCC state yields a 3-body 
interaction like in (fT6)l . while a Kitaev's code produces 
the standard 2-body classical Ising model, namely, 



Rising (/? J) 



(id) ' 



(23) 



ii/ Symmetry: A distinctive feature of our result ((T6|) is 
that the classical spin model associated to the TCC state 
does not posses the up-down Z2 spin-reversal symmetry. 
However, the partition function pB]) exhibits a Z2 x Z2 
symmetry. Recall that the lattice A is 3-colorable at sites, 
so that we can redundantly label our classical spin vari- 
ables as erf with c = r,g,b the color at site i. Then the 
change of variables 

-> s{c) at, s{r)s{g)s{h) = 1, s{c) = ±1, (24) 

gives a global symmetry, which has symmetry group 
Z2 X Z2 because s{b) = s{r)s{g). The ground states 
have to display this symmetry, in fact. Consider states 
in which the values of the spin variables only depend 
on the color, that is, for which erf ~ fc with fc = ±1. 
Such states can be labelled with the tag {fr,fg,fb)- 
Then it is easy to check that for the ferromagnetic case 
J > the ground states are labelled with the posi- 
tive parity tags (+, +, +), (+, -, -), (-, +, -), (-, -, +) 
whereas in the antifcrromagnetic case J < 
they arc labelled with the negative parity tags 
(-,-,-),(-,+,+),(+,-,+),(+,+,-). Thus, each 
parity sector, or classical ground state of (fT5|) is fourth- 
fold degenerate. Notice that the 3-body Ising model in 
such 3-colorable lattices of triangles shows no frustration, 
as opposed to the standard Ising model (|23p in such lat- 
tices which is indeed frustrated. 

Remarkably, the gauge group underlying the topologi- 
cal order related to Hamiltonian ([3]) is also Z2 x Z2. The 
situation is the same also with toric codes, where the 
global symmetry of the classical system is Z2 and the 
gauge group for the toric code topological order is Z2. 
This is certainly not a matter of chance, since one can 
relate the types of domain walls in the classical system 
to the types of condensed strings in the quantum system, 
iii/ Selfduality: The models in the triangular and 
Union Jack lattices turn out to be self-dual like the usual 
2-body Ising model, with a critical temperature Pc given 
by the same condition. 



sinh 2A' 



1, K. 



PcJc = 0.4407. 



(25) 



model like (fT6l) or (|23p . A high-temperature expansion 
is a polynomial in the variable u = tanh(/3J) that is 
small when T ^ 00, while a low temperature expan- 
sion is another polynomial in the variable u* := e"^'^"' 
that is small in the limit T — > 0. Then, a self-duality is 
a relationship between the high-tcmperature expansion 
of one classical spin model in a given lattice A and the 
low-temperature expansion of the same lattice. This is 



precisely the case of the 3-body Ising model ( 16 ) on both 

2I and the 



Duality is a property between high-temperature and 
low-temperature expansions of a statistical mechanical 



the triangular and Union Jack lattices [2l| 
standard Ising model ^ 

iv/ Universality Classes: Interestingly enough, the 
3-l3ody Ising model on the triangular lattice [l3|, [iBl 
and the Union- Jack lattice [l^ are exactly solvable mod- 
els under certain circumstances and this fact allows us 
to check their criticality properties when compared with 
those of the standard Ising model solution. 

The critical exponent for the specific heat in the 3-body 
Ising model on the triangular lattice is a = |. This rep- 
resents a power law behaviour which is in sharp contrast 
with the well-known logarithmic divergence {a = 0) of 
the specific heat in the standard Ising model (|23p . Other 
representative exponents are also different: the correla- 
tion length exponent is = | (vs. v = 1), the magneti- 
zation exponent is (3 ~ (vs. /9 = |), while they share 
the same two-point correlation function exponent at the 
critical point ?/ = -j. 

For the 3-body Ising model on the Union Jack lattice, 
the specific heat critical exponent is also remarkably dif- 
ferent a = ^. In fact, if the coupling constant J is allowed 
to be anisotropic, then even the critical exponent a may 
take on a set of continuous values in (0, ^) depending on 
a parameter related to the coupling constants [l9| . 

The computational capabilities of a topological color 
code depends on the 2-colex lattices where it is defined. 
For a TCC on a square-octogonal lattice it is possible 
to implement the whole Clifford group of unitary gates 
generated by the set of gates {H^K'i ,K2], where H is 
the Hadamard gate, AT 2 the 7r/8-gate and A2 the CNOT 
gate Q. However, for a 2-colex like the hexagonal lat- 
tice the set of available gates is more reduced since the 
7r/8-gate cannot be implemented topologically Q. Thus 
we point out a remarkable connection between different 
computational capabilities of color codes that correspond 
to qualitatively different unversality classes of their asso- 
ciated classical spin models, despite the fact that both 
color codes have the same topological order. 



C. Borders 

If we want to consider classical systems of spins with 
a finite number of sites, then we have to introduce either 
borders or a nontrivial topology. Since in TCCs the non- 
trivial topology gives rise to degeneracy, it is preferable 
to have borders. Also, borders play a role for the ideas 
to be explained in the next section. 
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Figure 6: Here both an hexagonal 2-col6x £2 and its dual 
triangular lattice A (dashed) are displayed to illustrate how 
borders are introduced in the 2-colex if A has borders. All 
vertices of the 2-colex which are not triangles in A have been 
removed, and also all the faces which keep no vertices. The 
faces that only keep part of their vertices remain, but only 
their Z face operators are kept in the stabilizer. 



In TCCs, borders can be of several types. For exam- 
ple, in it was shown how to build borders of a given 
color. Here our guide to construct the border must be 
the dual lattice A, which now has a border, along with 
the properties of the classical system. Then, as shown 
in Fig. [6l in order to construct the stabilizer for a TCC 
with border in such a way that (|17p remains true is to 
start with an infinite 2-colex £2 and then keep only part 
of it following certain criteria, (i) Keep those vertices 
of £2 which correspond to triangles in A. (ii) For those 
faces f of £2 which keep all their vertices, we keep the 
face operators X\ and Zj. (iii) For those faces f which 
only keep a subset f of their vertices, we introduce a 
face operator acting on those qubits. Condition (i) 
ensures the correspondence between triangle chains in A 
and string- nets in £2- Conditions (ii) and (iii) ensure the 
correspondence between Fq and Aq. 



IV. MEASUREMENT-BASED QUANTUM 
COMPUTATION WITH COLOR CODES 

In a measurement-based quantum computer (MCQ) 
0, information processing is carried out via a sequence 
of one-qubit measurements on an initialized entangled 
quantum register. This is an alternative to the stan- 
dard gate-based quantum computation that can simulate 
quantum networks efficiently. 

An interesting problem is to study the performance of 
the MQC when the initial entangled multiparticlc state 
is a topological toric code like in Kitacv's model 0]. In 
particular, under which general circumstances the MQC 



based on the planar Kitaev code can be efficiently simu- 
lated by a classical computer. The answer to this ques- 
tion is that the planar code state can be efficiently sim- 
ulated on a classical computer if at each step of MQC, 
the sets of measured and unmeasured qubits correspond 
to simply connected subsets of the lattice . 

Likewise, another very interesting problem is whether 
a topological color code state is classically sumulable in 
an scenario of measurement-based quantum computation 
(MQC). By extending the results of sect.lIIH it is possible 
to address this problem here. Our aim is to see what 
conclusions can we learn from the statistical mechanical 
models in order to test the classical simulability by MQC 
of the topological code states. 

The results of sect. IIIII can be interpreted as a com- 
plete projective measurement of the planar topological 
color code state (fTTjl onto a very specific product of sin- 
gle qubit states p2p . This type of global measurements 
are not enough for doing a MQC based on the color code 
state. Instead, we need to allow for more general one- 
qubit measurements and to perform them in an adaptive 
fashion as the computation proceeds from the starting 
point till the end. 

To be more specific, let us consider a generic qubit 
state with complex coefficients 

:=c2|0)„+cj|l)„, cS,cJ eC. (26) 

Then, a MQC starts with a planar color code (fTTjl and 
we apply a scries of projective measurements := 
\x)„{x\, a; = 0, 1, from the first qubit = 1 in the code un- 
til the last one = The order in which this sequence 
of measurements is carried out through the 2-colex lattice 
is arbitrary. After each measurement, the corresponding 
qubit at the vertex gets projected onto one the states 
in ([26]) and the result is a value for the outcome denoted 
by m„ = 0, 1. 

The result of a run of a MQC is a set of outputs 
nil, ■ • ■ , JTio, . . . , m\<x}\ with a certain probability distribu- 
tion P(mi, . . . , m|(jj|). We adopt the definition Q that a 
MQC is classically simulable in an efficient way if there 
exists a classical randomized algorithm that allows one 
to sample the outputs mi, . . . , m|(jj| from the probability 
distribution P(mi, . . . ,TO|(jj|) in a time poly(|2J|). 

Then, after a complete run of a MQC with pT|) we get 
the following generalized overlapping 

Omqc (*c|(g)|</5)«. (27) 

Using the same type of computations that led to (|14p . we 
arrive at the following expression 

Omqc = E l"^)" 

-n^sE n fe)- 

DG2J 7Gro DeaT:a;„=l ^ " 
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This result can also be turned into the partition function 
of a statistical model with 3-body Ising interactions (fT6|) 
but with complex and inhomogeneous Boltzmann weights 



Wijk 



(29) 



depending on each triangular plaquette < i,j,k >. The 
generalized overlapping ([^5]) is proportional to a gener- 
alized partition function of a 3-body Ising model with 
inhomogeneous and complex coupling constants 



Z(/3,{J,,fe}) := ^e''5:(.„.)eA-'..--^-.--, 



(30) 



where the lattice A can be the complete triangular or the 
complete Union Jack lattice. 

At an intermediate stage of a run of a MQC with the 
color code state, the 2-colex will split into two subsets 
£2 9Jl U DJl, with DJl the set of measured qubits and 
Tl the set of unmeasured qubits Thus, during the 
running of the MQC starting with the color code state, 
we shall find generalized partition functions of the type 
(|30|) but for a lattice that is the dual of the subset of 
measured qubits: A = 9Jl*. 

Therefore, we arrive at the situation that in order to 
classically simulate a topological color code state in a 
MQC we need to simulate the conditional probabilities 
P(mB+i|TOi, . . . ,TO|ti|) (at step 0-1-1 knowing the prob- 
abilities of previous steps) and for these we need to be 
able to compute efficiently in a time poly(L) on the size 
L of the intermediate lattice at step 0-1-1. 

At this point there is a sharp difference between the 
classical simulation with MQC of Kitaev states and color 
code states. Kitaev states can be classically simulated 
under very general conditions: the subsets 9Jl and dJl 
must be simply connected. The basic ingredient to 
achieve this result in the 2-body Ising model is that even 
though a generalized standard Ising model (with arbi- 
trary complex and inhomogeneous couplings) may not be 
translationally invariant, nevertheless there always exist 
a technique allowing it to be mapped onto a dimer cover- 
ing problem (DCP) which in turn can be solved efScie ntly 
through the Pafhan method in polynomial time 0, [23l |. 

However, the dimer problem technique is applicable to 
2-body interactions but not for the 3-body interactions 
that arise in the generalized statistical mechanical models 
from color codes. In the case of the 3-body Ising model 
with uniform and real couplings Jjji; = J e R in a tri- 
angular lattice, it can be exactly solved by mapping it 
onto the generating function of a suitable site-colouring 
problem (SCP) on a hexagonal lattice 17|, [l^, which 



can be solved by the Bethe ansatz. In order for this site- 
coloring mapping to work, certain restrictive conditions 
on the triangular lattice must be fulfilled. In particular, 
the partition function (jl6p has to be defined on a periodic 
triangular lattice with L rows in the horizontal direction 
and N columns in the vertical direction. Let us denote 
it as Z'^^jl^. Let us also denote by Zff^ the generating 



function of a site-coloring problem on a hexagonal lattice 
with M = |L and N columns. Then, the aforementioned 
mapping works in the limit N 



7(3) 



ySCP 



N-»l. 



(31) 



Furthermore, the SCP is solved by Bethe ansatz tech- 
nique. This technique also poses another fundamental 
problem in this situation since it is used to compute the 
eigenvalues of the associated transfer matrix of the SCP 
and then the issue about the completeness of that spec- 
trum in terms of Bethe ansatz eigenfunctions arises. This 
issue is always a difficult question and, strictly speaking, 
it is a conjeture. Quite on the contrary, these difhcul- 
ties are absent in the standard Ising model case since the 
DCP is more versatile. 

The situation becomes even more difficult if we con- 
sider the generalized partition function pop in the frame- 
work of an intermediate step in the MQC. Then, it is not 
known how to solve it efficently with a mapping to a SCP 
in an hexagonal lattice without restrictions. 

This site-coloring mapping plays a similar role than 
the dimer covering mapping in the standard 2-body Ising 
model. However, the known solutions to this site-coloring 
problem demand more restrictive conditions on the type 
of lattices and they are less powerful than the dimer map- 
ping technique. 

As for the topological color code on a 2-colex like the 
Union Jack lattice, similar conclusions apply: in the the 
case of real couplings Jijk = J G R, it is exactly solvable 
since it can be mapped onto a 8- vertex model [l9| , which 
in turn has to be solved by the Bethe ansatz. 

The fact that the 3-body classical Ising model is ex- 
actly solvable in the triangular and Union Jack lat- 
tices for real and isotropic couplings, is not enough so 
as to conclude that the corresponding topological color 
code states are classically simulable in a MQC scenario. 
Therefore, the classical simulability of topological color 
codes with MQC remains an open problem. 



V. CLUSTER STATES AND MODELS WITH 
AN EXTERNAL FIELD 

In a MQC scheme of quantum computation, the in- 
put state is called a cluster state Q which is a rather 
general entangled state associated to a great variety of 
graph states, i.e., states constructed from qubit states 
located at the vertices of a lattice specified by the inci- 
dence matrix of a graph. A very important property of 
these cluster states is that they can be created efficiently 
in any system with a quantum Ising-type interaction be- 
tween two-state particles in the specified lattice configu- 
ration. In Sect II VI we have assumed that the input state 
for the MQC was a TCC state, without caring about how 
it could be prepared. Here we show how such a topolog- 
ical color code state can be obtained from a appropriate 
cluster state. As a by-product, this construction will turn 
out to be useful for obtaining classical Ising models in 
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Figure 7: The graph needed to obtain a color code state from 
a cluster state. The graph is bipartite, and the vertices are 
divided in black and white. Black vertices correspond to the 
faces of the 2-colex. White vertices correspond to the vertices 
of the 2-colex. Black thick lines represent the edges of the 
graph, and grey lines correspond to the edges of the 2-colex. 

an external magnetic field being associated to color code 
states. 



u G 1X2 we will obtain a series of binary values = x^. 
The remaining qubits are then in a state characterized 
by the conditions: 

Vfe;?, ZflAi) =Xf|«'c(x)). (34) 

Thus, if Xf = for all the faces f, the result is the TCC 
state |\l/c) ([TT|) . In the other cases the result is essentially 
a TCC, in particular the state is 

|vl/,(x)) := J2 l7) (35) 

where x denotes the binary vector of the measurement 
results and Fx is the set of string-nets 7 with dcj = 

In fact, the original cluster state can be written as 
follows: 

where |x) = (^u^y, |a:^u)u is an element of the computa- 
tional basis of the subsystem of qubits in il2 . To check 
that this is indeed the correct expression for the cluster 
state, it is enough to note that it satisfies (|32p. 



A. Cluster state formulation of TCC 

Instead of giving a general definition of cluster states, 
we will consider only bipartite cluster states. So let & be 
a finite bipartite graph, that is, a graph such that its set 
of vertices il is the disjoint union of two sets il = ili Uil2 
in such a way that neighboring vertices never belong to 
the same ili. Consider the quantum system obtained by 
attaching a qubit to each of the vertices u. For each such 
vertex, we denote by N{u) the set containing both u and 
its neighbors. The cluster state \k) for the graph is 
then completely characterized by the conditions: 

Vueili, Xw(u)|k> = 

Vueil2, Zjv(u)|«;) = (32) 

In order to relate the TCC C of a 2-colex £2 to a cluster 
state, we construct a bipartite graph by setting ili = QJ 
and il2 = "S- Then the edges are defined so that Ui = G 
ill is a neighbor of U2 = f G il2 if D is a vertex of f in £2 , 
see Fig. [71 Observe that the corresponding cluster state 
\k) satisfies 

VfGj, X^\k) = \k,), (33) 

because = ^^^^X^f^p). In fact, if 7 is a string-net 
with (?c7 = 0, then = However, to keep things 

as simple as possible, let us just consider 2-colexes in 
which all closed string-nets are boundaries. If we measure 
in the Z basis all the qubits corresponding to vertices 



B. Models with an external field 

Thus far, we have only considered statistical mechan- 
ical models with zero external magnetic field. Here we 
go beyond that situation, considering the formulation of 
models with 3-body Ising interactions and arbitrary mag- 
netic fields from the projection of topological color codes 
onto appropriate product states. 

Let us define the product state 

|$p(J,h)} := (g) |0(J„))„ (g) |(/)(/if))j, (37) 
oeUi fGU2 

where J = ( J„) G RI*^'! , h = (h^) G RI^! and 

|0(s)) :=cosh.s|0)+sinhs|l), s G R. (38) 

Consider the overlapping 

0(/3,J,h):=(M'c|<i>p(/3J,/3h)). (39) 

Its value is 

5^(x| (g) i0(/.f))f E (ti (g) = 

l[cosh{(3h^)l[cosh{(3J,)J2 E nV'Tl""''"' 

f X jer^ f 

(40) 
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where 

:= tanh(/?/if), u„ ta.nh{f3J„). (41) 

We want to relate (p9)) to the partition function of 
a classical spin system. As in section IIIH we consider 
the lattice A dual to the 2-colex £2 and we associate a 
classical system to A by attaching classical spin variables 
CTj = ±1 to each of its sites i. This time however we 
include triangle dependent couplings Jijk in triangles and 
a site dependent external field hi. Thus, we want to 
derive the high temperature expansion for the partition 
function 

Z(AJ,h):=^e-^«(-''»'), (42) 
{"} 

where the classical Hamiltonian is 

H(J, h) - ^ /liO-j - ^ J^jkUtCTjCTk- (43) 

We start using the identities 

^ cosh(/3/i,) + a, sinh(/3/i,;), 
^l3J,,.a,a,a, ^ cosh(/3Jyfe) + a,a,ak sinh(/3J,,fc), (44) 

so that, 

1 + u^a.,) n (1 + UijkO-iCr-jak), 

(45) 

where 

C :=[|cosh(/3/i,) J] cosh(/3Jyfc), (46) 

:= tanh(/3ft,,), u,jk :^ tanh{pjijk). (47) 
Let us rewrite ([45|) in the form 

Z(/3,J,h)=C5^^5^n("^'^'')"" n i^^Jk<J^'JJ<TkY'^'' , 
X 5 {a} i {i,j,k) 

(48) 

where x = (x;) is a binary vector and S = 
S(iifc) ^ijk^ijk is a chain of triangles, that is, a formal 
sum over triangles with binary coefficients. Reordering 
the expression we get 

Z(/3,J,h) = c5^^e(x,<5)n«."' n ^^J-fc"""' 

X (5 i {i.j,k) 

(49) 

where 

6(x, -5) := E n n (^^'^.^fc)'-^" ' (50) 

{a} i {i,j,k) 

From (HIl) it follows that 

e(x,<5) = (f' !f^;na.).=-«:> (51) 
1^ 0, m other case. 



where (j, k)i denotes the pairs (j, k) which form a triangle 
with i, and iV is the number of sites. Finally we can give 
the desired high temperature expansion of the partition 
function: 

Z(/3,J,h)=2^c5] H"'"' n (52) 

X SeAx i (i.jM) 

where Ax contains those chains of triangles such that at 
any given site i an even (odd) number of triangles meet 
if Xi = {xi = 1). 

In order to compare (j40|) and (|52|) . we first observe 
that sites i correspond to faces of the 2-colex f S = il2 
and triangles Ay^ correspond to vertices of the 2-colex 
S QJ = ill. This correspondence relates in an obvi- 
ous way Xi with Xo, hi with h„ and so on and so forth. 
Also, there is an exact correspondence between chains of 
triangle 6 and string-nets 7. In particular, this correspon- 
dence identifies Ax with Fx, so that we get the desired 
relationship between the overlapping and the partition 
function 

2(/3,J,h) = 2^0(/3,J,h). (53) 



VI. CONCLUSIONS 

We have shown that the classical spin models asso- 
ciated to quantum topological color code states in the 
two dimensional lattices called 2-colexes are Ising models 
with 3-body interactions. We have studied this mapping 
in the triangular and the Union Jack lattices, which are 
the duals of the two very represenative 2-colexes, namely, 
the hexagonal and the square-octogonal lattices, respec- 
tively. This is a genuine difference with respect to the 
case with toric code states which yield the partition func- 
tion of the standard Ising model with two-body inter- 
actions. Ising models with different n-body interactions 
have very different properties in general. Remarkbly, dif- 
ferent computational capabilities of the topological color 
codes depending on the chosen 2-colex correspond to dif- 
ferent universality classes of the associated classical 3- 
body Ising models. 

The tools employed to relate classical spin models with 
topological color code states can be extended to study 
the performance of such topological states when they 
are considered as input in a measurement-based quan- 
tum computation. Then, the classical 3-body models 
that arise involve arbitrary complex couplings and lat- 
tice shapes. The problem of evaluating their correspond- 
ing generalized partition functions cannot be performed 
with the dimer covering technique that is so successful 
in the case of the classical 2-body Ising model. The 
similar technique in the 3-body case is a particular site- 
coloring problem that only in very specific instances has 
been solved by means of the Bethe ansatz. The com- 
pleteness of the Bethe ansatz poses in turn more fun- 
damental problems in this regard. Therefore, the fact 
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that the 3-body Ising model is exactly solvable in certain 
conditions is not enough to conclude so far that the par- 
ent quantum topological color color states are efBciently 
classically simulated by MQC. 

Another interesting result is the construction of a clus- 
ter state from which we can construct the topological 
color code state. This turns out to be useful in order to 
obtain classical spin models in the presence of arbitrary 
external magnetic fields. Likewise, there are other two- 



dimensional multipartite states that arise in the study of 
quantum antiferroma gne ts that may lead to a variety of 
classical spin models [25| . 
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